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Zero temperature crossover from superfluidity to superconductivity with carrier density increasing is studied for 
a two-dimensional system in the s-wave and d-wave pairing channels. It was assumed that the particle attraction 
correlation length depends on carrier density n/ as ro ~ 1/ ^/nj. Such a dependence was found experimentally for 
the radius of magnetic correlations in La2-xSrxCuOi. The short range Coulomb repulsion was also taken into 
account. It is shown that the behavior of the system with doping is fundamentally different from the case with 
ro{nf) = const. In particular, similarly to the d-wave case, the crossover in the s-channel takes place only if the 
coupling is larger of some minimal value, otherwise the Cooper pairing scenario takes place at any small carrier 
density. The relevance of the model to the high-temperature superconductors is discussed. 
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1. INTRODUCTION 

The problem of the crossover from superfluidity 
to superconductivity with charge carrier density 
or coupling constant changing has a long history 
jl]. The interest to this phenomena has arisen 
again after the discovery of high-temperature su- 
perconductors (HTSCs) in 1986 It was 
already known upon that time that the superflu- 
idity of composite bosons transforms into the su- 
perconductivity of overlapped Cooper pairs with 
chemical doping in the case of the s-wave pair- 
ing. Now, in the s-wave pairing case the problem 
is quite well explored for the 3D systems 
and, particularly for the quasi-2D case Q. For 
the 2D case this problem was studied at T = 
( when a long-range superconducting is still pos- 
sible in a 2D system [^) for the case of local 
attraction (see, for example 1^,^,^) and for the 
phonon-exchange model . Most of these prob- 
lems are reviewed in jl^ . 

Due to a layered structure and the anisotropic 
symmetry of the order parameter in HTSCs, the 
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crossover in the 2D system in the d-wave pair- 
ing channel is of a special interest. However, be- 
cause if its complexity, this case is not so well 
understood at present. The d-wave pairing for 
the case of the extended Hubbard model with 
the nearest neighbor attraction case was analyzed 
in [^|j9|. The crossover from superfluidity to su- 
perconductivity in the s-wave and c?-wave pairing 
channels for a 2D continuum model was studied 
in paper [|l0|, where also thermodynamic proper- 
ties of the system in the crossover particle density 
region were considered. The authors proposed an 
interaction potential, which is attractive at dis- 
tances between particles shorter of some value rg 
and longer of ri, and repulsive due to electron 
correlations at short distances, r < ri. It was 
found, in particular, that there exists a minimal 
value of the attractive coupling constant in the 
d-wave pairing channel, which gives the crossover 
from superfluidity to superconductivity at small 
carrier densities, i.e. the fermion chemical poten- 
tial changes its sign and becomes positive with 
carrier density growth. 

The correlation length tq was assumed in 10 
to be a parameter, which does not depend on the 
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carrier density Uf per cell. However, concern- 
ing HTSCs it would be interesting to consider a 
more realistic case, when at small carrier densities 
ro ~ l/y^/nj and the proportionality coefficient is 
of order of a lattice constant. Such a dependence 
was experimentally observed in HTSC cuprates 
for the length of spin-spin correlations, which are 
believed to be responsible for the hole attraction 
in these materials. In particular, it was found for 
La2-xSrxCuOi in underdoped regime, that the 
magnetic correlation length decreases with car- 
rier density per cell according to the dependence 
3.8A/^ 

In what follows we analyze the possibility of low 
carrier density crossover for a model analogous to 
with particle repulsion at distances r < ri 
and attraction at ri < r < ro, where, however. 



ten in a standard form 



(1) 



(a is parameter of order of the lattice constant, 
and its possible value is discussed in the next Sec- 
tion). Obviously, at large carrier densities such 
that ri > ro the superconductivity in this sys- 
tem should disappear. However, since we are in- 
terested in the small carrier concentrations, it is 
assumed that this relation does not take place. 

As it will be shown below such a dependence 
of ro leads to a qualitatively different behav- 
ior of the system with respect to the case with 
ro(n/) = const. In particular, there exists a mini- 
mal value for the coupling constant when the two- 
particle bound states exist at low carrier densities 
in the s-wave pairing channel. The existence of 
such a threshold value of the coupling constant is 
typical for the d-wave pairing case (see, for exam- 
ple 1 10 ij]). Another interesting property is: for 
any coupling constant in both channels there ex- 
ists a corresponding carrier density value, below 
which the system is in the superconducting state. 



2. THE MODEL AND THE MAIN 
EQUATIONS 

The Hamiltonian of the system which describes 
the non-retarded fermion interaction can be writ- 



dxidx2V{xi,X2)ip]{xi)ipl{x2)ipi{x2)->jj]{xi), (2 



where m is the effective fermion mass, and fi is the 
chemical potential; fermi-operators ip,y{x) depend 
on the space-time coordinate x — (x, t) . The in- 
stantaneous interaction potential is chosen in the 
next form 



V{xi,X2)^6{h-t2)V{r), 



with 



V{r) = VrepOin - r) - VattMr - n)9{ro - r), (3) 

which corresponds to potential used in [ p^ . Here 
r = |xi — X2I is an inter-particle distance. Posi- 
tive parameters Vrep and Vattr correspond to par- 
ticle repulsion at r < ri and particle attraction 
at ri < r < ro. The charge carrier density depen- 
dence of the correlation radius ro is given by (]l|) 
with the parameter a = A/2/7rao, ao is the square 
lattice constant. This relation can be easily es- 
timated from the equality (7r/2)rgAf/ — aQNceii, 
where on the left side the volume of the 2D sys- 
tem is expressed as a volume (circle of the ra- 
dius ~ ro) occupied by one particle, multiplied 
by the full number of particles TV/, iV'^'^" is an 
elementary cell number in the system. The free 
fermion bandwidth W is connected with ao as 
W = 7r^/(mag). It should be noted, that the re- 
lation (|^) at a = ^J^JiiaQ is in a good agreement 
with the experimental data for La2-xSrxCuO\ 
p3[ | , where the plane magnetically ordered lattice 
parameters are equal to 5. 354A and 5.40lA, and 
the corresponding parameter a is ~ 3.8A. 

In order to study the superconducting proper- 
ties of the model in the channels with different 
pair angular momentum Z, it is convenient, sim- 
ilarly to [pO| , to approximate the Fourier trans- 
form of (tJl) by a separable potential: 



k, k. 



-A/W;(ki)w; (k2), 



(4) 



where A; is an effective coupling constant, and 
functions 



w;(k) = hi{k)cos{lipk) 



(5) 
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with coefficients 



(6) 



(i + fc//coy+i/2' 

k is the momentum modulus k = |k|, and 
(pk is its angle in polar coordinates k = 
{kcos{(py:), ksin{(p]s)). Parameters ko and ki put 
the momentum range in the proper region. They 
are connected with the potential (|^) parameters 
as fco ~ and fci 1/ri. Below we put 

fco = 1/fo and ki ~ 1/ri. Obviously, the ex- 
pression is not the exact Fourier transform of 
but it sets the interaction in right momen- 
tum range and has the correct asymptotic behav- 
ior at small and large momenta: V^^^^^ ~ ^1^2 
and ^jk2 ~ 1/Vfci^2, correspondingly. Since we 
are mainly interested in the low carrier density 
region, where the crossover can take place, the 
correct behavior of the interaction potential at 
small momenta is the most important. In the 
case of low carrier concentrations the small mo- 
menta give the main contribution in the integral 
of the gap equation (see Eq. (0) below). We shall 
study s- and d-wave channels with I = and 2 
separately, so we assume that the parameters A; 
for both channels are independent. It will be as- 
sumed also that the pairing takes place at zero 
pair momentum at T = 0. 

The minimization of the ground state energy 
with respect to the superconducting order param- 
eter A/(a;i - X2) = V{xi - X2) < ^■i{xi)'4>iix2) > 
and the chemical potential gives in the case of the 
approximation (Q) the standard pair of equation 
for the crossover problem (see, for example |p^ ): 

dp A,(p) 



A,(k) = -A, J 
dp 



{2ny 
where e(p) 



(1- 



(27r)2 2V£2(p) + A,(p)2 



(7) 



Ve2(p) + A,(p)^ 



(8) 



p^/2m — /i is a free particle disper- 



sion law. The form of the equation (|^) allows to 
search the solution for the superconducting order 
parameter in the next form 

Ai(k) = AoiU.z(k), (9) 

where Aq; does not depend on momentum k and 
wi (k) is defined by (^) and (||) . In this case the so- 
lution of the system (0),(^ gives the dependence 
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Figure 1 . a) Crossover line gs — sp is presented 
for the s-pairing channel (solid line). The dotted 
line represents the corresponding curve for the 
case with r-o(n/) = const at ro — uq. Here and 
below all parameters are expressed in units of the 
bandwidth W. b) The same as a) at very low 
charge carrier densities. 
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of the gap parameter Aq; (Aqs for I = and 
Aod for I = 2) and the corresponding chemical 
potential fi on the particle density n/ and the 
coupling A/. 

3. CROSSOVER FROM SUPERFLUID- 
ITY TO SUPERCONDUCTIVITY 

3.1. s-wave pairing channel 

The solution of the system (0),(||) at jj. = 
gives the crossover line As(n/), which separates 
the parameters regions, where the local (/i < 0) 
and Cooper pairing > 0) take place. Note, 
that in the s-wave pairing case the Coulomb re- 
pulsion parameter ki does not enter in the equa- 
tions. The numerical solution for the crossover 
line for the s-wave pairing case is presented in 
Fig.l, where we have put Fermi energy ep instead 
of Uf, since in the 2D case they are connected by 
a linear relation ep = f / {ma^) , and the dimen- 
sionless coupling constant gs = mXs/{iiT) is intro- 
duced. The solution for the ro{nf) = const case 
is also presented. As it follows from the numerical 
calculations, there exists a minimal value of cou- 
pling =~ 2.574, which corresponds to carrier 
density e^? ~ 6.8 x 10~''VF on the crossover line, 
necessary to generate the crossover from super- 
fluidity to superconductivity with doping, other- 
wise the Cooper pairing regime takes place at any 
charge carrier density. In other words, there is a 
minimal value of coupling constant which leads 
to the two-particle bound states in the s-wave 
pairing channel at small ep- It is important to 
note, that there is no such a minimal coupling 
in the s-wave pairing channel when the correla- 
tion length ro(n/) = const Moreover, the 
crossover with charge carrier density increasing 
in the s-wave channel takes place for any known 
doping independent interaction potential |l^ ]. 

Another interesting property, which follows 
from the Fig. lb), is the "inverse" crossover from 
superfluidity to superconductivity with charge 
carrier density decreasing at small values of ep 
{sp < 6.8 X 10~^VK). This is a consequence of 
the competition between two opposite processes 
which occur with ep decreasing. The lowering of 
Sp tends the system to become a superfluid, but 
at the same time it leads to growing of and 
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Figure 2. The dependence of Aqs on ej;- at dif- 
ferent coupling parameters gs is presented for the 
s-wave case. 
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Figure 3. The dependence of ^ onsp is presented 
for the s-wave pairing channel at different values 
of gs. 
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coupling constant at low doping 
9i 



I 1 21 1+1/2 



(10) 
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Figure 4. The crossover line g^isp) is presented 
for the d-wave case (solid line). The dotted line 
is the crossover curve for the case ro(ri/) = const 
at ro — uq. 



makes the pair size larger, i.e. the pairs become 
bounded weaker. On the other hand, the carrier 
localization on dopants at very small carrier den- 
sities also makes a possible density region of su- 
perfluidity more narrow |T^. However, the region 
of extremely small particle densities sp is 
not very interesting from the point of view of the 
connection with HTSC, since in this region the re- 
lation ro ~ 1/ y/nj does not hold, and the model 
is not correct. As it will be shown in the next 
Subsection, in the d-wave case this crossover from 
from superfluidity to superconductivity with dop- 
ing decreasing takes place at rather large Fermi 
energy values. 

The doping dependencies of the gap Aqs and 
the chemical potential n at different values of the 
coupling constant gs are presented in Figs. 2, 
and 3. The gap is increasing with the doping 
almost linearly, except the region of extremely 
low concentrations. It is interesting to note, that 
the chemical potential in the superfluidity region 
(/i < 0) has a minimum at a finite value of ep and 
it is equal to zero at e^^ = 0, since the effective 



is zero at £i? = 0. In other words, the local 
pairs are the most strongly coupled at some fi- 
nite carrier density value. This is also a conse- 
quence of the competition between the Fermi sur- 
face formation and correlation radius decreasing 
with charge carrier density growth. This situa- 
tion is qualitatively different from the ro{nf) = 
const case, where fco does not depend on parti- 
cle density and renormalized coupling constant 
is £i?-independent in the s-channel. In this case 
fj, — Ei,/2 < at e_F = 0, where i^f, < is a 
two-particle bound state energy. 

3.2. d-wave pairing channel 

In this case the Coulomb repulsion parameter 
fci is present in the equations (0) and How- 
ever, the presence of this parameter leads just 
to renormalization of the dimensionless coupling 
constant g^ = mA(i/(47r) and the energy gap pa- 
rameter: 55 = gd/iei/Wy, Ag^ = Aod/ei, where 
£1 = k\/{2m) is characteristic energy of Coulomb 
repulsion. We shall consider the case, when the 
Coulomb repulsion is much smaller than W , i.e. 
this is the case of a large free-fermion bandwidth. 
Therefore, the coupling constant is assumed to 
be large: gl - 10'' - 10^ (see Figs. 4-6). The 
crossover line g^iiep) for the d-wave pairing case 
is presented in Fig. 4. Qualitatively, the behavior 
of the system with doping and coupling changing 
in the d-wave pairing channel is similar to the s- 
wave case. The important difference is that the 
low carrier density superconducting state exists 
at rather high values of e^. Also in this case 
there is a minimal value of the coupling constant 
for two-particle bound states g^ ~ 1.7856 x lO'' at 
Ei? ~ 0.2731VF. It should be noted that the exis- 
tence of the large threshold value for the coupling 
constant in both channels can be a possible an- 
swer on the question why the crossover has been 
not observed in cuprates. 

The charge carrier density dependence of Aq^ 
and fi at different coupling parameters g^ are 
presented in Figs. 5 and 6. The superconduct- 
ing gap Aod at low charge carrier densities is 
much smaller than in the s-case, due to stronger 
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eF-dependence of effective coupling constant at 
low carrier densities (|l^) . The magnitude of the 
order parameter Ag^ starts to grow almost lin- 
early with sp increasing when the Fermi energy 
is larger of some minimal value (Fig. 5). This be- 
havior is qualitatively similar to the doping de- 
pendence of the gap of cuprates in the under- 
doped regime. However, because of its simplicity, 
the model can't describe the decreasing of the gap 
with charge carrier density increasing at large val- 
ues of HTSCs. For this other properties of the 
charge carrier interaction in cuprates have to be 
taken into account (see the last Section). 

The small-£ f superconducting region with /i > 
as well as region of superfluidity are rather large 
in the d-wave pairing channel. Also in this case 
the chemical potential is equal to zero at £f = 0, 
i.e. there are no two-particle bound states at very 
low charge carrier densities. It should be noted 
that the decreasing of the chemical potential with 
increasing charge carrier density at small e p both 
in s- and d-wave pairing cases indicates a nega- 
tive electronic compressibility. This can be re- 
lated to increasing of antiferromagnetic correla- 
tions at low carrier densities. 

4. CONCLUSION 

Theoretical description of the behavior of su- 
perconductor with carrier density changing is an 
interesting and important problem, in particular, 
because of its possible association with HTSCs. 
The microscopic mechanism of the superconduc- 
tivity in cuprates is not known so far, and the 
solution of phenomenological models, which take 

into account some of the properties of HTSCs, 
can help to clarify the nature of their unusual 
behavior, and maybe even help to understand 
the microscopic mechanism of the HTSC phe- 
nomenon. 

In this paper the possibility of the crossover 
from superfluidity to superconductivity with 
charge carrier density and coupling constant 
changing in different pairing channels at T = 
was studied for a model, which qualitatively takes 
into account one of the properties of HTSCs, 
namely, the doping dependence of correlation 
length ro at low carrier densities. It has been 
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Figure 5. The dependence of AJ,^ on ep is pre- 
sented for the d-wave pairing channel at different 
values of coupling constant g^. 
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Figure 6. Chemical potential fj. as a function of 
Ei^ is presented for the d-wave case at different 
values of coupling constant g^. 
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shown, that even this simple model results in 
interesting and unusual properties, which are 
rather different from a more standard case with 
ro(n/) = const. In particular, the two-particle 
bound states in the s-channel exist only if the 
coupling constant is larger of the threshold value, 
similar to the d-wave pairing case. At any value 
of coupling constant larger of the threshold one, 
the "inverse" crossover from superfluidity to su- 
perconductivity takes place with doping decreas- 
ing in both s- and rf-wave channels. Of course, 
such a simple model can not pretend to describe 
doping dependence of the gap and chemical po- 
tential of HTSCs. The momentum dependence of 
the interaction potential has to be taken into ac- 
count more carefully, especially in the overdoped 
regime, where £f is rather large and the separa- 
ble potential may be not correct. In general, also 
the effect of the retardation of interaction can not 
be neglected. These and some other questions are 
planed to be studied in a future work. 
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